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ON THE ZEROS OF THE PEARCEY INTEGRAL AND A 
RAYLEIGH-TYPE EQUATION 


GERARDO HERNANDEZ-DEL-VALLE 


Abstract. In this work we find a sequence of functions /„ at which 
the integral 

(1) v(t,x) = r 

J —oo 

is identically zero for all f > 0, that is 

v(f,/„(f))=0 Vt>0. 

The function v, after proper change of variables and rotation of the path 
of integration, is known as the Pearcey integral or Pearcey function, in¬ 
distinctly. We also show that each f„ is expressed in terms of a sec¬ 
ond order non-linear ODE, which turns out to be of the Rayleigh-type. 
Eurthermore the initial conditions, which uniquely determine each /„, 
depend on the zeros of an Airy function of order 4 defined as 

As a byproduct of these facts, we develop a methodology to find a 
class of functions which solve the moving boundary problem of the heat 
equation. To this end, we make use of generalized Airy functions, which 
in some particular cases fall within the category of functions with infin¬ 
itely many real zeros, studied by Polya. 
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1. Introduction 

The Pearcey integral was first evaluated numerically by Pearcey ifTTI in 
his investigation of the electromagnetic field near a cusp. The integral ap¬ 
pears also in optics flU, in the asymptotics of special functions [|6l, in prob¬ 
ability theory IfTSlI . as the generating function of heat (and hence Hermite) 
polynomials of order Ak for k gN ITTBI . It also falls into the category of 
functions considered by Polya IfT^ . that is functions with countably many 
zeros. For the numerical evaluation of the zeros of the Pearcey integral see 
for instance O, this will be important since this zeros will correspond to 
the initial value of each function /„(0). 

The main motivation in finding the zeros of the Pearcey function, which 
solves the heat equation 

( 2 ) ht{t,x) = 

is due to the fact that the main building block used to construct the density 
of the first time that a Wiener process hits a boundary /, is to find a function 
/ such that 

h{t,f{t))=0 Vt>0. 

For example, suppose there is a financial contract which will be activated 
if ever the price of an asset S (modelled as Brownian motion) reaches a 
prescribed boundary /. For instance, in Figure the blue line represents 
the evolution of the price of St, for t G [0,10], in turn the red line represents 
a boundary which activates a contract if it is ever reached. In particular, the 
barrier option is a contract of this type. For a more detailed exposition see 
for instance (Si . 

Next, we note that for some constant G M, the function 

with {t,x) G X M, solves the heat equation ([^. This is true since h in 
@ is a linear combination of the fundamental solution of @ and its first 
derivative with respect to the space variable x. It is clear that the function 
Q equals zero at x = —bt. Hence, for any a G M, and setting x = a —bt in 
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Figure 1 . The dotted line represents the priee of some as¬ 
set S, in turn, the solid line represents a boundary whieh will 
aetivate a eontraet the first time it is reaehed. In this exam¬ 
ple the eontraet was aetivated approximately at time t = 5. 
This is a random time, sinee it would have been impossible 
to foresee the outeome. 


Q we obtain 



We note that the right-hand side of this identity is in faet the density of the 
first time that a Brownian motion hits a linear boundary [|71 p. 196]. In 
praetiee these results are used for instanee in (a) the valuation problem of 
finaneial assets, in partieular in the valuation of barrier options [see Bjork 
(2009)], (b) in the quantifieaetion of eounterparty risk [see Davis and Pisto- 
rious (2010)], and in general in physieal problems. 

The main eontributions of this work are, on the one hand, finding the ze¬ 
ros of the Peareey integral. On the other, advaneing in the direetion of de¬ 
veloping a rather simple and straightforward methodology to find explicit 
solutions of the time-varying boundary problem for the heat equation. In 
this regard we note that there exist techniques to study the latter aforemen¬ 
tioned problem in terms of solutions to integral equations [|4]|. We recall 
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that solutions in terms of integral equations, in general ean only be evalu¬ 
ated numerieally. In turn, our approaeh leads to solutions in terms of ODEs. 

The paper is organized as follows. In Seetionj^we introduee the Airy 
funetion of order 4. Next in Seetion we define the Peareey integral and 
deseribe its eonneetion with the Airy funetion of order 4. In Seetion]^ we 
derive a Rayleigh-type equation, whose solution kills the Peareey funetion. 
The teehniques deseribed in Seetion]^ are illustrated with examples in See¬ 
tion 1^ In Seetion we derive a funetion whieh asymptotieally solves the 
moving boundary problem for the Peareey integral. This approximation ean 
be helpful in the numerieal solution of the Rayleigh equation. We eonelude 
in Seetionwith some final remarks. 


2. Generalized Airy function of order 4 

With respeet to the zeros of Fourier integrals, Polya proved ifT^ that all 
the zeros of 

/ oo _ 

g-M form = 1,2,3... 

-oo 

are real and infinitely many for m > 1. In turn, the generalized Airy funetion 
^ of order 4 ean be expressed as a solution of the following ODE 

(5) 0(3) = Y0, 

== (y — -hY0(-^~3) for 7 > 3. 

One ean prove, for instanee applying the Fourier transform to Q and solv¬ 
ing the resulting equation, that 0 is a partieular ease of Q when m = 2, 
namely 



Furthermore funetion 0 is symmetrie, with eountably many zeros in the real 
line—henee oseillatory—and tends to zero as it inereases to ±oo, see Figure 
1^ Regarding the zeros of Q, there exist asymptotie estimates whieh are 
derived by means of the method of steepest deseent [|T4ll . 
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Figure 2. We plot function (p, defined in Q. We observe 
the function is even and its first zero is at ±2.44197. 

3. The Pearcey integral 
Paris dH [IS analyzed the asymptotic behavior of 

/ oo 

Jm, n g N, r > 2, 

-oo 

which by rotation of the path of integration [u = te^) and use of Jordan’s 
lemma (see [fT4l f can be expressed as 

P'A.y) = P,{x.y) = eS r 

J — OO 

with X = and y = Ye^. In particular, the Pearcey integral, which 

solves Q, is the case P 2 . More explicitly, we have the following. 

Definition 3.1. WH The Pearcey integral is defined as 

/ CX3 

-00 

In this work we study instead the following Fourier integral 



(8) 
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because the zeros of v, for {t^x) G M+ X M, are expressed in terms of a 
continuously differentiable function /, as opposed to Q. See [fT4ll . 

Remark 3.2. We observe that the function v in ([^ is the convolution be¬ 
tween the kernel of standard Brownian motion and the generalized Airy 
function of order 4, in equation ©• 


4. Zeros of the Pearcey function 

Remark 4.1. Throughout this work, the n-th partial differentiation with 
respect to the space variable x of any given function v{t,x) is denoted as 
vW. 

In this section we find the function / for which the Pearcey function is 
zero for every t >0. The idea is to exploit, on the one hand, the differential 
form of the Airy function of order 4, defined in Q, and on the other to use 
the fact that the Pearcey function solves the heat equation Q. 

The main result is the following. 

Theorem 4.2. Suppose that v is as in ©, ^ solves ©, ^ is such that ^(^) = 
0, and f is a solution to the following Rayleigh-type ODE 

(9) /"(<) = 2 

with /(O) = ^ and f{ff) — /[2(j)^^'> {^)]. Then, for every t>Q,we 

have 

v{tj{t))=0. 

Proof Given that (j) is as in Q, its Fourier transform equals 

^(A) = exp 

Furthermore, if we apply the Fourier transform directly to the ODE Q we 
have that 
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Since this expression is already in Fourier domain, we eonvolve the previous 
expression with the heat kernel as follows 

This in turn, and by direet applieation of the integration by parts formula, 
yields 

(10) = xv + tv^^\ 

as well as 

(11) = V+Yv(^)+tv(2), 

after differentiation with respeet to x. 

Next given that there exists an /, see Polya IfT^ . sueh that the following 
holds for all t 

v{tj{t))=0, 

we differentiate (Leibniz integral rule) v with x = f{t), defined in ([^, with 
respeet to t to obtain 

1 poo C 1 

2i ~ --Xh--\dX = 0, 

whieh is equivalent to 

(12) /'(»)»'<"(»,/(<))+ ivP>(,,/(»)) = 0 

and 

(13) /"Wv(i)+/(t)(/(t)v(2)+v(3)) + iv(4) = 0 

after differentiation with respeet to t twiee. We note that equations pO] ) and 
( fTTj ) obtained from the Airy differential equation ([^ , as well as ( [T^ and 
( [T3| ) obtained from the heat equation, involve derivatives of v up to order 
4. What remains is to obtain Q from these expressions. To this end, from 
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( fT^ and ( [TT] ) we first have 

vi'Ht.m 

»■<'"(»,/(<)) _ , »■«(»./(<)) 
v(‘>(«,/W) ''^' v(l)(,./(,))' 

Next from ( [T^ and ( [T3] ) it follows that 


/"(<)+/'(») 


/'(‘)+ 


2 v(»(,,/(,)) 


0 


[•'' v('>(t,/(t))J 4 /(<))■ 


These identities yield 

/"(<)+/'(<) (- 2 (/'(<))" + <)=-^(/(')- 2 </'(<)). 

This completes the proof of Theorem |4.2[ □ 


5. Examples 


To illustrate Theorem 4.2 we next present examples and numerical ex¬ 
periments. 


Numerical Example 5.1. One can show that at t = 0 the following two 
identities hold for f: /(O) = ^ = 2.44197 and f\0) = 0.729925. Hence, 
from Theorem ^~2\ we may plot the solution of^, see Figure^for t G [0,4]. 
The ODE was solved numerically using Mathematica. 


Airy4[x_] : = 

1/(2* Pi)* NInte grate [Exp [I*x*y — y''4/4], (y,—Infinity, Infinity}] 
Airy4I [x_] : = 

NIntegrate[I Exp [ I x y — y''4/4]*y, {y, —Infinity , Infinity )]/( 

2 \[Pi]) 

Airy42[x_] : = 

NIntegrate [ —Exp [ I x y — y''4/4] * y''2 , {y, —Infinity , Infinity )]/( 

2 \[Pi]) 

xl = FindRoot [Re [ Airy4 [x ] ] , (x, 2.44)] 
xl = x /. xl 

fp = -N[Re[Airy42[N[xl ]]]]/(2*Re[Airy41 [N[xl ]]]) 
s = NDSolve [ ( g ’ ’ [ x] == -(g[x]/4) - 
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Figure 3. Numerical solution of 0 with /(O) = ^ = 
2.44197, and f'{0) = 0.729925. The graph was plotted us¬ 
ing Mathematica. 

1/2 X g’[x] + 2 (g’[x])''3, g[0] == xl, 
g’[0] == fp } , g, {x, 0, 4}, AccuracyGoal —> 20, 

PrecisionGoal —> 10, WorkingPrecision —> 33] 

Plot [ Evaluate [g [x] /. s], {x, 0, 4}, PlotRange —> All] 


Numerical Example 5.2. To test the accuracy of the solution in Example 


t e s t [ x_ ] : = g [ X ] /. s 

F0[t_, x_] : = 

1/(2* Pi)* NIntegrate [ 

Exp[I*x*y — y''2*t/2 — y''4/4], {y, —Infinity , Infinity)] 
Table[Re[FO[i/100, t e s t [ i / 1 0 0 ][ [ 1 ] ] ] ] , {i, 0, 400, 1}] 

Next, we present some further examples of the methodology discussed in 
the previous section. 

Example 5.3. The Airy function of order 3 solves the following ODE 

= x(j){x). 

Using the same argument as in the proof of Theorem \4.2\ it follows that 
v^^\t,x) = xv{t,x) +tv^^\t,x). 


5.1 we may use the following code in Mathematica in the interval t G [0,4]. 
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This last expression evaluated at f yields 


v(2) 

7I) <'•/(»)) = »■ 


Finally, from the previous expression and ( [7^ it follows that 

, 1 v(2) 


2 ' 


Thus, for some constant C 


See [[m p.l26 ]. For applications of the Airy function in the first hitting 
time problem of Brownian motion up to a quadratic function see IlSlI . 


Example 5.4. Given the following ODE 

(l)^'^\x) =x(j){x) +(l)^^\x), 

and following the same line of reasoning as in the proof of Theorem \4.2\ we 
have 


v^'^\t,x) = xv{t,x) -\-tv^^\t,x) + 


From {12) we obtain 


It follows that 


v(2) 


fit) 


t 1 


or 

t t^ 

/W=-2-4+C 

For example if C = —2.58811 the moving boundary problem of the heat 
equation associated with ^ is solved. 
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Example 5.5. Given the following Bessel ODE 



Similar calculations as in the previous examples yield 


fit) 


Example 5.6. The derivative of the Airy function Aifx) solves 


X(j)"ix) = (p'ix) +x^(j)ix), 


which yields 


—tv^^\t,x) + (t^ —x)v^^\t,x) + (1 +2xt)v^^\t^x) + + i ) v (?, y ) = 0 . 


Alternatively, from Example \5.3\ we also have that 

v^^\t,x) —tv^'^\t,x) — a) — 2v{t,x) = 0. 


Using the same arguments as those described in Section^leads to 


2f{t)f'{t) + {l+tf{t))=0. 


(14) 


This is the Abel equation of the second kind and its solution can be ex¬ 
pressed in terms of the Airy function of order 3 Ai, as follows: 



See Eigure^for a numerical example with /(O) = —1.01879. 


6. Zeros of the Pearcey function. An asymptotic approach 

In this section we carry out an analysis in order to find an asymptotic solu¬ 
tion to the moving boundary problem associated with the Pearcey integral. 
This result is useful when finding the numerical solution of the Rayleigh 
equation Q. The main result of this section is the following. 

Theorem 6.1. Suppose that v is as in (|^, ^ is any zero of the Airy function 
of order 3, and 


/(t) = -2(t/3)3/2 + ^(3t)i/6, t>0. 


(15) 
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Figure 4. The Numerical solution of /, defined in Exam¬ 
ple 5.6 and such that /(O) = —1.01879 was plotted with 
Mathematica. 


Then 


as t ^ oo. 




Proof. For brevity let us just consider the term within the brackets in 
i.e., 

expj/Ay-^A^t-^j. 


Introduce a variable 


A4 


exp |A(/y-|- at) — -f a)^t “ ^ 

Set M = A -f a and rearrange terms to obtain 

exp ^{iy+at) — \ my 

f 1 

= exp < M(/y-|- at-l- a^) — -M^(t-l-3a^) -l- au^ — — 


/Y^ 1 p . 

^-\-^ah~iay _ 
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To get rid of the heat (or quadratie) term note that 


( 16 ) 

That is, 


t + = 0 , 


gives 


a± = ±i\j 


exp I ui 




33/2 


,3/2 


± 




31/2 


Next, if we ehoose a+, as in ( [T^ , 


exp |mz 




33/2 


A/2 


31/2 “ 4 ' 


and thus, from @, 

^ ^ 2kJ-oo 

Now, let z = u{3t)^/^, u = z(3t)^^/^, and {3t)^^/^dz — du, whieh yields 


(30 

or equivalently 




1 

2;r, 




(30^6 ^+'"3 4i3,)V^dz, 


/ 5 2 /r I r 

{3ty/^e-T6‘^-^3yv{t,y) = ^ e ( 30 '/^ ^ dz. 

27r J —00 


Letting y = —2{t/3)^/'^-\-^{3tY/^ 








Henee for arbitrary ^ 

f 1 //C >3 £ f2/3 / r#-|3/2 

In partieular if is a zero of the Airy funetion then 


+ ^(30’ 


1/6 I 


lm|(30^/'’e^"^?^v(^L-2 


3/2 


+ < 3(30 


1/6 


= Ai(^). 

= Ai(^) 
= 0 . 


□ 
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Numerical Example 6.2. For instance, if ^ = —2.33811 we have that 
f{t)=2{t/3)^/^ + 2.33Sn{3ty/^. SeeFigure^ 

Next, using the results of this Seetion and Seetion]^ we present an algo¬ 
rithm whieh ean be used to solve @ in Theorem |4.2| . 

Algorithm 6.3. For some root ^ of the Airy function Ai, let f be defined as 
in (0- Given an arbitrary time t G M’*' we may find a solution to equation 
0 in the interval [t — efi -\- e] as follows. 

Given t and find a root of function v, defined in (|^, starting at f{t). 
Lef fit) be the root obtained in the previous step, 
then fit) = -\/2-v^^\tJit))/v^^\tJit)). 

With fit) and fft) solve ^ in [t — £,t + e], for some £ > 0. 

Numerical Example 6.4. Suppose we choose t = 10 and £ = 2. Then the 
procedure is the following 

Fl[t_, x_] : = 

1/(2* Pi NIntegrate [(I*y)* 

Exp[l*x*y — y^2/2*t — y^4/4], {y, —Infinity , Infinity}] 

F2[t_, x_] : = 

1/(2* Pi J* NIntegrate I( I*y }'^2* 

Exp[I*x*y — y'^2/2*t — y^4/4], (y, —Infinity , Infinity}] 
asy[t_] := 2*(t/3)^(3/2) - AiryAiZero[I]*(3* t)^(1/6) 
xIO = FindRoot[Re[F0[10, t]], (t, Nfasy[10]]}] 
xIO = t /. xIO 

fpIO = -Re]F2]I0, xI0]]/(2*Re]FI]I0, xIO]]) 

sIO = NDSolve ] {g ’ ' ] X ] == —l/4*g]x] + 2*(g’]x])^3 — l/2*x*g’]x], 

g]I0] == xIO, g’]lO] == fplO}, g, (x, 8, 12}, AccuracyGoal -> 20, 
PrecisionGoal —> 10, Working Precision —> 301 

1 . Possible applications and work in progress 

Due to the stochastic and periodic nature of several economic variables, 
as for instance Mexico’s general CPI or the Fruit and Vegetable annual in¬ 
flation and assuming VF is a random walk, these processes can be modelled 
as 

(17) fij smi(j)j + 2nvjt) + t = 1,2,..., 

;=i 






ON THE ZEROS OE THE PEARCEY INTEGRAL 


15 



Figure 5. The red line is the real boundary /, the blue is 
as in with ^ — —2.33811. The graph was plotted with 
Mathematiea. 

where the /3/ and represent respeetively the amplitude and phase at a time 
given frequeney Vj. In turn a eontinuous time approximation of ( [T7| ) ean be 
expressed in terms of the solution of an SDE of the form sinee 

(18) dXt = ^{t,Xt)dt + o{t,Xt)dWt, 


AXj = ^U.Xj)Aj + aU.Xj)AWj, 7 = 0 , 1 ,..., 

where functions /i and o are respectively: 

n 

^ 2;rv//3jCos(<))j + 2;rv/), and o{t,Xt) = l, 

7=1 

and VF is a Wiener process. A reasonable set of questions that could be 
asked could be for instance: 

What is the probability that Fruit and Vegetable annual 
inflation will reach 20 points before the end of 2016? 

What is probability that the general CPI will remain 
between 3 and 4 percent until the end of 2017? 
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As it turns out, to answer the previous questions it is necessary to under¬ 
stand the moving boundary problem of heat equation addressed in this work. 
More specific examples is still work in progress. 

8. Concluding remarks 

In this work we find the zeros of the Pearcey function, in terms of the 
solution of a Rayleigh-type equation. This goal is achieved by exploiting, 
on the one hand, the differential equation of an Airy function of order 4 
and on the other by using the fact that the Pearcey function is a solution 
of the heat equation. As a by-product we develop a methodology, using 
straightforward techniques, to solve the moving boundary problem of the 
heat equation in the case in which the convolving function is a generalized 
Airy function. We expect that the techniques described within can be used 
in the construction of densities of the first hitting time problem of Brownian 
motion. The scope and applicability to the latter problem is still work in 
progress. 
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